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We apply the Ermakov-Lewis procedure to the one-parameter damped modes y recently introduced by Rosu and Reyes, 
which are related to the common Newtonian free damping modes y by the general Riccati solution [H.C. Rosu and M. Reyes, 
Phys. Rev. E 57, 4850 (1998)]. In particular, we calculate and plot the angle quantities of this approach that can help to 
distinguish these modes from the common y modes. 



PACS number(s): 03.20.+i 

In a previous paper hereafter denoted as I [Q , the non- 
uniqueness of the factorization of linear second-order dif- 
ferential operators has been exploited on the example of 
the classical Newtonian free damped oscillator, i.e. 



Ny 



dt 2 



2/3 



dt 



y = Q 



(1) 



The coefficient 2/3 is the friction constant per unit mass 
and o>o is the natural frequency of the oscillator (SI units 
assumed all over the work). 

The more general supersymmetric partner equation 
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d?_ 
dt 2 



2(3 



27 2 



dt 



y = o (2) 



1 ( 7 *+l) 2 , 

has been obtained in I. This new second-order linear 
damping equation contains the additional last term with 
respect to its initial partner (1), which may be thought of 
as the general Darboux transform part of the frequency 
jlj . T = I/7 occurs as a new time scale in the Newtonian 
damping problem. If this time scale is infinite, the ordi- 
nary free damping is recovered unless for the critical case 
which is special even in ordinary damping. As explained 
in I, the y modes can be obtained from the y modes by 
operatorial means. In the following we shall call them 7 
modes. For the three types of free damping they have 
been obtained in I as follows: 



(i) For underdamping, /3 2 < 



\J to 2 , — P 2 the underdamped 7 modes are 



denoting u v 
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lo u sva.{uj u t 



7< + 1 



cos(u; n i + (f>) 



(ii) For overdamping, 1 > J 2 and uj = ^ffi 2 



(3) 
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the overdamped 7 modes are 



yo = —A e 131 uj smh(io t 



(iii) For critical damping, /3 2 — 
solutions are given by 



cosh(o; i + <f>) 

(4) 

The critical 7 



Vc 



-An 
7* + l 



1 



-0t 



(5) 



These are the only possible types of one-parameter 
damping modes related to the free damping ones by 
means of Witten's supersymmetric scheme Q and the 
general Riccati solution ||. 

In practice the new parameter 7 can be very close to 
zero. In this case, it is very difficult to differentiate the 
7 modes from the ordinary ones. The only means we 
can think of is by recording somehow the geometric an- 
gle associated to the 7 modes and compare it with the 
same quantity in the ordinary damping cases. One is led 
to this conclusion noticing that the 7 modes have time- 
dependent frequencies uj 2 (t) = uj 2 — ^^^2 an d hence 
for them the Ermakov-Lewis (EL) procedure can be natu- 
rally applied || (for a recent review, see ||). For iv ^ /3, 
Eq. (2) can be reduced to a Bessel equation and the so- 
lutions can be written as follows 



and 



= r 1 / a 



AJ3 (fcr) + BY3 (kr) 



-1/2 



CIs(kT) + DKs.(kr) 



-0T 



(6) 



(7) 



where t = -ft + 1 and k 2 = 



When k — > 00 (i.e., 



7 — > 0), we can do Hankel's asymptotic expansions, i.e., 
of large Bessel argument but fixed Bessel order (we shall 
not reproduce these formulas here, the reader is directed 
to Abramowitz and Stegun 0). The point is that one is 
indeed able to get the solutions obtained by operatorial 
means from inspecting Hankel's expansions. Thus, the 
supersymmetric operatorial procedure gives merely the 
asymptotic 7 — > solutions, which however could be the 
most relevant from the physical viewpoint in this context. 

In the EL approach the angular quantities are given 
by the following formulas PJlO[] 



A0 d 



-2/3*' 



e 2 ^ p 2 



dt (8) 



and 



A6 S 



r(e 2 ^ pp) 



dt 
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(9) 
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for the dynamical and geometrical angles, respectively. 
Thus, the total angle will be 



fT p-2/3* 

A0 l = I -^dt 



(10) 

The so-called Pinney function p is the solution of Pin- 
ney's nonlinear equation || 



p"(t)+p(t)p'(t)+q(t)p = 



C 

pHt) 



exp 



-2 / p(t')dt' 



(11) 



2^ 



for C = constant (=1), p{t) = 2(3 and q{t) = lJq — jz^iy 
For p ^ constant there is a definite prescription of cal- 
culating p in terms of two independent solutions of the 
corresponding linear equation. We have followed the 
method of Eliezer and Gray || for p(t) in terms of lin- 
ear combinations of the aforementioned Bessel functions 
(for A = B = C = D = 1) that satisfy the initial condi- 
tions as given by those authors. In the critical damping 
case, we used the modes of Eq. (5) with A c = D c = 1. 
The results of the calculations for some particular val- 
ues of the parameters are plotted in Figs. la,b,c, 2a,b,c, 
3a,b,c for the 7 underdamped, overdamped, and critical 
cases, respectively. For comparison, the angle quantities 
for 7 = 0, within the same calculational scheme, arc dis- 
played in Figs. la',b',c', 2a',b',c', 3a',b',c', respectfully. 
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Fig. la 

The dynamical angle in the underdamped case for the 
following set of parameters: ujq = V%, f3 = 1, 7 = 0.1. 
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Fig. lb 

The geometric angle in the underdamped case and the same 
parameters. 




Fig. lc 

The total angle in the underdamped case and the same 
parameters. 
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Fig. 2a 

The dynamical angle in the overdamped case for loq — 1, 
f3 = y/2, 7 = 0.1. 



Fig. 2c 

The total angle in the overdamping case for the same 
parameters. 




Fig. 2b 

The geometric angle in the overdamping case for the same 
parameters. 



Fig. 3a 

The dynamical angle in the critical case for u)o 
7 = 0.1. 



= /3 = 1 and 




too, P parameters as in Fig. la and 7 = 0. 



Fig. 3b 

The geometrical angle in the critical case for the same 
parameters. 




Fig. 3c 

The total angle in the critical case for the same parameters. 



Fig. lb' 

The geometrical angle in the underdamped case for the same 
luo, f3 parameters and 7 = 0. 
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Fig. lc' 

The total angle in the underdamped case for the same luo, 
parameters and 7 = 0. 
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Fig. 2a' 

The dynamical angle in the overdamped case for the same 
loo, (3 parameters as in Fig. 2a and 7 = 0. 



Fig. 2c' 

The total angle in the overdamped case for the same luq, P 
parameters and 7 = 0. 




Fig. 2b' 

The geometrical angle in the overdamped case for the same 
too, P parameters and 7 = 0. 




Fig. 3a' 

The dynamical angle in the critical case for luq = f3 = 1 and 
7 = 0. 
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Fig. 3b' 

The geometrical angle in the critical case for loq 

7 = 0. 



= (3 = 1 




Fig. 3c' 

The total angle in the critical case for luq — (3 = 1 and 7 = 



